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Abstract 

We show the existence of thermal rectification in the graded mass quantum chain of har- 
monic oscillators with self-consistent reservoirs. Our analytical study allows us to identify 
the ingredients leading to the effect. The presence of rectification in this effective, simple 
model (representing graded mass materials, systems that may be constructed in practice) 
indicates that rectification in graded mass quantum systems may be an ubiquitous phe- 
nomenon. Moreover, as the classical version of this model does not present rectification, our 
results show that, here, rectification is a direct result of the quantum statistics. 
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A fundamental challenge in statistical physics is the derivation of macroscopic phenomeno- 
logical laws of thermodynamic transport from the underlying microscopic Hamiltonian systems. 
However, after decades, a first-principle derivation of the Fourier's law of heat conduction, for 
instance, is still a puzzle [I]. Many works have been devoted to the theme [2], most of them by 
means of computer simulations. But, in these problems where the central question involves the 
convergence or divergence of thermal conductivity, sometimes there is a significant difficulty to 
arrive at precise conclusions from numerical results [3j. Thus, the necessity of analytical studies 
together with the huge complexity of the associated nonlinear dynamical systems led to several 
works considering approximative schemes or simplified models [4j - since the pioneering work of 
Debye, the microscopic models to describe heat conduction are mainly given by systems of an- 
harmonic oscillators, leading to problems without analytical solutions. Many works involve, e.g., 
the use of approximations such as Boltzmann equations, master equations for effective models, 
the analysis of Green-Kubo formula, etc. An example of simple (effective) model that can be 
analytically studied is the classical or quantum harmonic chain of oscillators with self-consistent 
reservoirs. It has been proposed a while ago [5j [6], but it is always revisited [8|. In such 
model, each site of the chain is coupled to a reservoir; the first and last sites are coupled to 
"real" thermal baths, while the inner reservoirs only mimic the absent anharmonic interactions. 
The self-consistent condition means that there is no heat flow between an inner reservoir and 
its site in the steady state: the inner baths act only as a phonon scattering mechanism, such as 
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the on-site anharmonic potential in more elaborated systems. This model, the classical and also 
the quantum version, in opposition to a standard harmonic system, obeys the Fourier's law [7J. 

In this scenario of intensive study of the heat mechanism, problems involving the possibility 
to control the heat flow by constructing thermal nano-instruments such as thermal diodes, and 
even transistors, thermal gates and memories, have recently attracted considerable theoretical 
[TUl [TT] and experimental |12| interest. A thermal diode, or rectifier, is a device in which the 
magnitude of the heat current changes as we invert the device between two thermal baths. That 
is, in a thermal diode, heat flows preferably in one direction. There are some analytical attempts 
to explain the phenomenon of thermal rectification and/or design a diode by using simple models 
(see e.g. the spin-boson nanojunction [13], [H], and the billiard system [15]) but, again, most of 
the works consider numerical computations [9j [TD1 [TT] . It is worth to recall the extensive work 
of B. Li and collaborators: e.g. thermal rectification in asymmetric graphene ribbons is studied 
by using molecular dynamics simulations in ref. [16]; also in carbon nanocone structures |17j ; 
at silicon-amorphous polyethylene interfaces |18j ; in carbon nanotube intramolecular junctions 
[19], etc. A commonly used design of diodes is given by the sequential coupling of two or three 
chains with different anharmonic potentials [9| 110} ITT]. It is frequently studied, but it is also 
criticized due to the difficulty to be constructed in practice [ID]. Recently, Chang et al. [12], 
considering a different procedure, built a nanoscale thermal rectifier in an experimental work: 
they use a graded material, namely, nanotubes externally and inhomogeneously mass-loaded 
with heavy molecules, a system with asymmetric- axial thermal conductance. Graded materials 
are also considered in other numerical studies of diodes |20[I21). and the interesting results made 
some authors suggest that thermal rectification is guaranteed by the use of graded materials 
with anharmonic interactions. 

Having in mind the investigation of this suggested relation between graded (anharmonic) 
materials and rectification, in a recent paper [22] we study the classic harmonic chain of oscil- 
lators with self-consistent reservoirs (CSC). We show the absence of thermal rectification for 
the CSC model with graded masses or graded interparticle interactions by using the method 
developed in refs.[23], [23]. In [25], Segal obtains similar results by using the approach of [8]. 
We, and other authors [U [7J, understand the SC models as effective anharmonic models (the 
inner reservoirs mimic the absent anharmonic potentials), but it is not a consensus. Anyway, 
even if we discard any relation with anharmonicity, these results show, at least, that the onset of 
Fourier's law does not guarantee thermal rectification in asymmetric classic models - we recall 
that the Fourier's law does not hold in standard harmonic systems (i.e., harmonic models with 
reservoirs at the boundaries only), but it holds in the SC models. 

Still searching for the mechanism behind the rectification in graded mass systems, we recall 
that, at low temperatures, quantum effects may introduce significant changes. The thermal 
conductivity for the quantum self-consistent harmonic chain (QSC), for example, depends on 
temperature [U [26], in opposition to the CSC thermal conductivity [7j. Thus, searching for 
possible quantum effects in the mechanism of thermal rectification, in a quite recent work [27], we 
turn to the inhomogeneous QCS and, in the linear response regime, i.e., for the chain submitted 
to a very small gradient of temperature and considering only linear corrections in computations 
with changes in the temperature, we show the absence of thermal rectification despite quantum 
effects in the conductivity. However, we may still ask about the possibility of a rectification 
for the quantum system submitted to a large gradient of temperature or, at least, for some 
effect beyond linear corrections. It is, indeed, a very important question: in ref. [281, the authors 
show that there is a significative rectification in some chaotic billiard systems "provided the 
temperatures (of the two sides of the system) are strongly different..." 
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In the present work, with the focus on the onset of rectification in inhomogeneous systems, 
we revisit the graded mass QSC model, and study it by using an analytical approach valid for 
any temperature gradient, i.e., beyond the linear response regime. For this quantum model, we 
show that, in opposition to the behavior of its classical version, there is a thermal rectification. 
We offer an explanation for this difference: in the quantum version, the expression for the heat 
flow ([3]) involves a distribution for the phonon frequencies (the frequencies in the system) that 
depends on temperature ©, leading to a mix of temperature, frequencies and also particle 
masses (see eqs.([3l S])). It does not happen in the classical system - see the comments for the 
high temperature limit at the end of this manuscript. This intricate mix of temperature, masses 
and frequencies makes the heat current change as we invert the graded mass system between 
two thermal reservoirs. The presence of rectification in this quite simple, say, bare model is 
of considerable interest: it indicates the generality of such phenomenon in the experimentally 
realizable class of graded mass systems, i.e., it indicates that rectification is not related to 
intricate interactions in the microscopic systems. Moreover, as far as we know, this work is the 
first demonstration of a "purely" quantum thermal rectifier, in the sense that, for such model, 
the rectification is a direct result of the quantum statistics: i.e., rectification is absent in the 
classical version of the model. 

Let us introduce the QSC model. We use a Ford-Kac-Mazur approach, as detailed presented 
in ref.[8], in order to describe the quantum system and its time evolution to the steady state. 
Here, all the baths connected to the chain are modeled as mechanical harmonic systems, with 
initial coordinates and momenta determined by some statistical distribution. Then, we solve 
the quantum dynamics given by Heisenberg equations, take the stochastic distribution for the 
initial coordinates of the baths, as well as the limit t — > oo, and obtain the expression for the 
heat flow in the steady state. 

The Hamiltonian of our system, a chain (W) with harmonic interparticle and on-site poten- 
tials, with each site connected to a bath (B), also with harmonic interactions, is given by 

n = u w + Y.' H B i + Y J x wVB i x Bi , (i) 

i=l i=l 

Us = \x T s M s Xs + \x^^ s Xs, 

where S = W or Bi] M\y, M Bi are the particle-mass diagonal matrices for the chain and 
baths; &w an d ar e symmetric matrices describing the interparticle and on-site harmonic 
interactions; and V Bi gives the interaction between the z-site and its bath (more details below). 
We have, for each part, X = [Xi, X2, • • • , Xn s ] t , where X r is the position operator of the r-th 
particle; X = M~ 1 P, where Pi is the momentum operator of the Z-th particle. Of course, it 
follows that [X r ,Pj] = ih5 rt i. Finally, the dynamics is given by the Heisenberg equations 

M W X W = -$ w X w -^2vlX Bi , (2) 

i 

M Bt X Bi = -* Bi X Bi -Y l VZ i X w . 

i 

The formulas for the heat currents inside the chain and from each reservoir to the chain 
are related to (XwX^) and (X B Xy^). Before presenting the formulas, let us give a very short 
resume of their derivation. To find the expressions, we turn to the Heisenberg equations (|2|), 
treat the equations of the baths as linear inhomogeneous equations, and plug these solutions 
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into the equations for the chain. Then, we take the average over the initial conditions of the 
baths, which are assumed to be distributed according to equilibrium phonon distributions with 
properly chosen temperatures, determined such that the self-consistent condition holds, that is, 
we must take the temperatures of the inner baths such that there is no heat flow between an 
inner bath and its site. We reach the steady state by taking the limit t — > oo. For technical 
reasons we still take to ~~ > —oo, and consider the Fourier transform of t. We note that, by 
plugging the solutions of the baths back into the equations of motion for the system, we get 
a quantum Langevin equation. Remember that, for our model, all particles are connected to 
heat reservoirs that are taken to be Ohmic. The coupling strength to the reservoirs is controlled 
by the dissipation constant £ defined by an expression involving the matrix V presented in the 
Hamiltonian above (see ref.|8j). 

The expression for the heat flow from the Z-th reservoir to the chain is given by 



N 



Ec 2 

771=1 



du to 



Lm 



(3) 



x — [/(o;,r,)-/(a;,r m )], 

7T 

where £ is the dissipation constant; 

[G w (u)}- 1 = -oj 2 M w + $w-Yl S / + ( w )> 



(4) 



the matrix X/" above has only one non-vanishing element: 
distribution for the Z-th bath 



H,l 



iQuj; f(u,Ti) is the phonon 



f(oj,T l ) = l/[exp(tvuj/KT l )-l}; (5) 

the variable u in the expressions above comes from the Fourier transform, namely, X(uj) = 
(27r) _1 f_°° dt X(t) e lwt . We stress that a detailed derivation of the heat flow formula (|3|) is 
presented in ref.jl]. 

The heat flow in the chain may be given by J-\, i.e., by the heat flow from the first reservoir, 
together with the self-consistent condition Ti = T$ = ■ ■ ■ = 7jv-i = 0, which means that the 
inner reservoirs do not inject energy into the system (the energy comes from the first reservoir, 
passes trough the chain and goes out to the last reservoir). In such condition, we have T\ = —J-n- 
To investigate the presence or absence of thermal rectification, we need to examine the heat 
current in a chain with reservoirs at temperatures T\,T2, ■ ■ ■ , Tjv, where T\ and TJv are given, 
and T2, ■ ■ ■ , Tjv-i are determined by the self-consistent condition [Ti = 0, . . . , -F/v-i = 0); and 
also in the same chain with inverted thermal reservoirs, precisely, with the baths taken at 
temperatures T[,T^ ■ ■ ■ , T' N , where T[ = Try, T' N = T\, and T^, . . . , T' N-1 are determined by the 
self-consistent condition (J^ = 0, . . . ,7-^^ = 0). In the case of absence of thermal rectification, 
we shall obtain T\ + J-[ = 0. That is what we will examine. 

The determination of Gw(w), the matrix given by eq. (jlj) and necessary to study Fi, may be 
a very difficult task - we know a precise solution for some specific cases: e.g., for a homogeneous 
next-neighbors interparticle interaction and all particles with the same mass M\y = ml [8] ; and 
also for the case of particles with alternate masses, i.e. rrij = mi for j odd, and mj = m2 for 
j even [26]. Hence, we first restrict our interparticle interaction $ to the homogeneous nearest- 
neighbor case, i.e., $ = —A, the lattice Laplacian, and study the equations for the smallest 
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possible chain: N = 3 (it is, in fact, a cell of a large chain). Then, we show that a graded mass 
distribution leads to a thermal rectification. After that, we argue to show that similar relations 
extend to N = 4, and so on. 
Thus, for A = 3, we have 



G9» 



Ri.i -l 
-1 R 2 ,2 -1 

-1 #3,3 



(6) 



where R kjk = 9 k k = % — w2 M k — iQw. To have graded masses, we take M\ < M 2 < M 3 . And 
using the notation c = ( 2 h/ir, A im = \(G w (w))i im \ 2 , fk = f(w,T k ), we get (from eq.©) 

T 3 

It is easy to see that Aj k = A k j. By taking the inverse of G^ above, we have 



c J dw w 3 A 12 [fr -f 2 ]+ cj dw w 3 A 13 [f 1 - / 3 ], 
c [ dw w 3 A 31 [f 3 -f 1 } + c fdw w 3 A 32 [f 3 - f 2 \. 



A 12 
A32 



-1|2 



[(2 -w 2 M 3 y + ( 2 w 2 ]/\detG wl , 
^-wH^f + CWyidetG^l 2 . 



(7) 



Note that, by using F\ = —F 3 , from the equations above we obtain (as expected) F 2 = 0. For 
the chain with inverted reservoirs, we obtain the expressions for F[ and F 3 similar to those for 
F\ and F 3 above, but with the change T\ f-> T3, i.e., with the replacement /1 0/3, and f 2 — > f 2 . 
The question is: do we have F\ + F[ = 0? 

From J 7 ! = —F 3 (i.e., = 0) and F[ = — F 3 (i.e., J 7 ^ = 0, we obtain 



J dw w 3 A 12 [fx -f 2 } = J dw w 3 A 32 [f 2 - f 3 ] , 

J dw w 3 A 12 [f 3 - f 2 \ = j dw w 3 A 32 [& ~ h}- 

To get the absence of rectification, we must have T\ + F[ = 0, i.e., 

J dw w 3 A l2 [h -f 2 ] = J dw w 3 A l2 [f 2 - h], 

or, which comes directly from eqs. (|8l9H0p (and also from F 3 + F 3 = 0) 

J dw w 3 A 32 [f 3 -f 2 } = J dw w 3 A 32 [f 2 ' - h\- 

That is, if we have eqs. (|8l9|10p . then eq. lfTTj) follows; similarly, given eqs. (|8|9llip . we have eq. lfTO 
Let us use the short notation 

Fi{T k ) -- 



(8) 
(9) 

(10) 
(11) 



dw w 6 A 12 f k , F 2 (T k ) = / dw w 6 A 32 f k 



and similarly for Fj(T 2 ). From eqs.([8|) and ([9]) we obtain 

: Fx(T{) + F 2 (T 3 ), 
■ F 2 {T 1 )+F 1 {T 3 ). 



Fi(T 2 ) + F 2 (T 2 ) 
F l (T^) + F 2 (T^) 



(12) 
(13) 
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We write the condition for absence of rectification (llOp as 



Fx(T 2 ) + iq(T 2 ') = ^(TO + F X (T 3 ), (14) 
or, which comes from these equations, 

F 2 (T 2 ) + T 2 (T^) = F 2 (Ti) + T 2 (T 3 ). (15) 

Let us briefly give an idea of our analysis before turning to the mathematical expressions. It 
follows that any function Fj(T) above, and also the sum Fi(T) + F 2 (T), is a monotone, strictly 
increasing, convex function of T (it is easily verified by studying its derivative d/dT). Hence, 
given Ti and T3, the value of T 2 is univocally determined by eq. (fl~2l) . and the value of T' 2 by 
eq. (|13p . Moreover, F±(T) is bigger and increases faster than T 2 (T), for M\ < M 2 < M3. In 
short, we have three different equations for T 2 and T 2 , that are determined by the two first 
equations; as the third one may not be derived (it is independent) from the first and second 
equations, we shall have rectification, i.e. the third equation will not be always satisfied. For 
linear equations in T (that is the case for the approximations of high temperature or linear 
response regime, as we show ahead), the system of two variables and three equations has a 
solution only if it is linearly dependent - more details ahead. 

Now we make explicit, i.e., we give the mathematical details of the analysis described above 
for the eqs. (|12l [T3l I14p . For simplicity, we will consider small temperatures (and so, also small 
changes in T), but, of course, we will go beyond linear corrections. Since F (F\ and F 2 ) is an 
analytic function of T (for T > 0), we write, for some Tq 7^ 0, \Tq — T\ small, 



F\(r) = F 1 (To) + F 1 / (ro)(r-r ) + ^-^(r-r ( 



0) 



Fl (T )-F[(T )T + ^l 



+ 



F[(T )-2T Q 



rp + F l( T 0) T 2 



= ao + aiT + a 2 T z , (16) 

that is, we discard corrections up to 0((T — To) 3 ). We have (at least for small To) ai,Q 2 > 0. 
The same follows for T 2 (T), 

F 2 (T) = b Q + b 1 T + b 2 T 2 . (17) 
Hence, introducing eqs. (fl~6l [T7]) in eq. (fT2l) . we get, for small T\ and T3, 

_ _ -(ai + 61) ± V(ai + ^i) 2 + 4(^i+5 3 )(a 2 + 6 2 ) 



2(a 2 + b 2 ) 

where 



(18) 



A 1 = aiTt + a 2 Tt, B 3 = hT 3 + b 2 T^ (19) 

Of course, only T 2 > makes sense, and so, T 2 is determined by the positive expression in 
eq. (fT8j) above. Similarly, for T 2 we get, from eq. (fT3|) . 



_ -(ai + 61) ± ^(ai + fri) 2 + 4(^1 3 + B 1 ){a 2 + 6,) 

1 2 — ttfz nn 



2(a 2 + 62) 

where 



(20) 



A 3 = aiT 3 + a 2 T£, B 1 = b 1 T 1 + b 2 Tl (21) 
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Again, we consider the positive expression for T 2 . 

Now we check if eq. (|14j) is verified or not, that means absence or existence of rectification 
respectively. Namely, we check if it is true the equality 

(a x T 2 + a 2 T 2 ) + {a x T' 2 + a 2 T'\) = {a x T x + atf?) + ( ai T 3 + a 2 T 3 2 ). (22) 

To make the investigation, we introduce the (positive) expressions for T 2 (|18p and T 2 ()20p in 
equation above. By taking T x = T, T3 = T + e, with small e, and using 

VI + x = 1 + -x - \x 2 + 0{x 3 ), 

we easily verify that eq. (|22p is not satisfied. E.g., by taking in eq. (|22p the terms involving only 
e and e 2 (not those with T, Te, etc), we get, for the L.H.S. and R.H.S. of eq . (|22l . 

L.H.S. = aie + -^V(°2 + & 2 )e 2 + (a 2 b x - a x b 2 ) ^ + ^ ~ 2 



a x + bi (ai + 61) 3 

R.H.S. = aie + a2e 2 . 

Note that there is an equality up to first order in e, and so, as we have already recalled [27], we 
do not see rectification in the linear regime. But, considering the corrections beyond the linear 
approximation, the difference is evident, which shows the existence of rectification. 

For transparency, let us also give more details of our analysis for some situations of absence 
of rectification. It happens, e.g., in the following cases. (I) For the system with homogeneous 
mass distribution (obviously): if M\ = M 2 = M3, then A X2 = A% 2 (and so, eq.fJHD is obtained 
by the sum of eqs. (|12p and (|13p ). (II) In the high temperature regime: in such case, as shown 
in ref.[8], the quantum system behaves as the classical one, where there is no rectification (as 
already shown in refs.[22j and |25j). To see this result from our previous equations, note that for 
large T we have w/[exp(hw/KT) — 1] — > KT/K, and so, eqs. ([T2"T) . (fT3|) . and (fb^) become linear 
equations in T 2 , T 2 , T\ and T3. It is easy to see that they are linearly dependent (and so, eq. (fT4|) 
comes from eqs. (H2p and (|13p ). (Ill) In the linear response regime (already shown in ref. |27| and 
commented above): again, as described in ref. [8], for small temperature gradients we have 



dw w 3 — (w,T)^2 A km{T} 



m=l 

where T above is some average temperature. The expression above leads us to the same case of 
(II) (i.e., we have linear equations in T 2 , T 2 , etc). 

Now we turn to larger chains. For iV = 4, we have the equations J-\ = — T± and T 2 = 0, 
^3 = (the self-consistent equation). But, as F\ + T 2 + T% + = (it may be easily verified 
from the expressions for J 7 ), if we take T\ = —T^ and T 2 = 0, then F% = follows (i.e., it is 
not a new equation; and vice- versa: T 2 = and J-3 = leads to F\ = — T£). Given T% and 
T4, these two equations (e.g., F\ = —Fa and T 2 = 0) determines T 2 and T3. The same follows 
for the system with inverted reservoirs (T 2 and T3). Hence, the equation for the absence of 
rectification T\ + J-[ = is an extra equation involving variables already determined, and the 
analysis follows as before (for N = 4 and any other N). We emphasize that the study of the 
dependence of the rectification on the system size (not presented here) is of importance: recall 
that, in the two-segment diode, the rectification vanishes as the system size goes to infinite. 
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We note that, for a system with three sites, if we take a mass distribution with M\ ^ M3 
(e.g., M\ = M2 7^ M3, or Mi 7^ M2 = M3) we still have rectification. For larger N, the analysis 
for the graded system is "automatic", but it may be intricate for any other mass distribution. 
It is an interesting (and difficult) problem to investigate which inhomogeneousity in the mass 
distribution is enough to guarantee the rectification. 

We still recall that, for a standard harmonic system, i.e., for the harmonic chain with reser- 
voirs at the boundaries only (in other words, turning off the inner reservoirs, say, the "schematic 
anharmonicity" ) , we have (see ref.[8j), 



and similarly for J 7 ', with the replacement f\ f-> /jv; i.e., T\ = —J~[, and so, there is no 
rectification. 

To conclude, our results show the existence of thermal rectification in the graded mass QSC, 
a quite simple model. It indicates that the presence of some phonon scattering mechanism, or 
some effective anharmonicity, is enough to guarantee the thermal rectification in an inhomoge- 
neous quantum model. We stress the importance of the quantum nature of the reservoirs that 
provides an energy mixing ([5]). It is probable that some other intricate mechanism, such as 
real anharmonic potentials, may provide such mix in classical models, as indicated by numer- 
ical studies [20]. Finally, we recall that graded mass models are not only theoretical systems: 
they may be constructed in practice [12]. In short, diodes of graded mass systems sound to be 
ubiquitous (and experimentally reliable) structures. 
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